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Abstract 

We compute the particle multiphcities and transverse energies at central and nearly 
central AA collisions at RHIC and LHC. The initial state is computed from pertur- 
bative QCD supplemented by the conjecture of saturation of produced partons. The 
expansion stage is described in terms of hydrodynamics assuming longitudinal boost 
invariance and azimuthal symmetry. Transverse flow effects, a realistic list of hadrons 
and resonance decays are included. Comparison with the data of the multiplicities at 
^/s = 56 AGeV and 130 AGeV from RHIC is done and predictions for the full RHIC 
energy and LHC energy are made for the multiplicities and transverse energies. The 
reduction from the initially released minijet transverse energy to the Et in the final 
state is less than in the one-dimensional case but still dramatic: a factor of 2.7 at 
RHIC, and 3.6 at the LHC. 
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1 Introduction 



The idea of parton saturation, the growth of the number of partons becoming inhibited 
when a partonic system reaches a sufficient density, was introduced some 20 years ago 
for pp colhsions . In the context of ultrarelativistic heavy ion colhsions, the idea of 
parton saturation in the wave functions of the colhding nuclei governing and regulating 
the final state particle production was discussed first in 0. It was suggested that all 
particle production in high energy AA collisions could be computed in perturbative 
QCD (pQCD) as the saturation scale Psat ~ ^^^^ becomes large, Psat ^ Aqcd- Also 
a pQCD approach supplemented by a soft QCD component |^ ^ implied that the 
perturbative mechanism would be clearly dominant in transverse energy production at 
high energies. Particle production in AA collisions has also been modeled in terms of 
classical gluon fields and in this approach the initial state parton saturation plays 
a key role as well. 

Saturation of produced partons at a perturbative scale psat can, however, be reached 
even if the partons in the wave functions of the colliding nuclei are not saturated. A 
consistent framework of pQCD parton production combined with the requirement of 
the saturation of the produced partons was introduced in 0, referred to as EKRT here. 
Based on the pQCD cross sections, it was shown that in central collisions the minijets 
with transverse momenta px > Psat, indicating an average transverse size of T^/plst 
for each minijet, are produced into the central rapidity unit so abundantly that the 
available transverse area, vri?^, is totally filled. Overcrowding thus takes place and the 
saturation of the produced partons occurs at the scale psat which grows with y/s and 
A. Most importantly, at RHIC and LHC energies and for large nuclei [A ~ 200), the 
saturation happens at perturbative scales, Psat ~ 1 ... 2 GeV ^ Aqcd- We will argue 
below that the average features of the produced minijet system are similar to those of 
a thermal system. Based on this, we will assume that the production time, determined 
by the saturation scale as 1/psat can also be taken as the formation time Tj of the 
(approximately) thermal QGP. With this assumption and the saturation hypothesis 
the initial state of the produced QGP can be calculated using pQCD alone. Especially 
in estimating the initial energy density, e = E^^^^ /{ixRlnAY), one can make use of 
the recent progress in computation of the minijet Ex production in next-to-leading 
order pQCD 0. 

At saturation the original minijet momentum distribution is already lost - as the 
produced gluons are interacting with each other - and the system has started to ther- 
malize. It was observed in that the produced, saturated, minijet system, which is 
~ 90% gluonic, looks indeed thermal from the point of view of the number of partons 
and the energy per particle. Therefore it is plausible to assume that the momentum 
distributions are close to thermal at saturation, and hydrodynamics applies in describ- 
ing the expansion stage of the system. Assuming an entropy conserving expansion 
stage, scaling laws with -^i and A were predicted in for the charged particle mul- 
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tiplicities. These predictions agree remarkably well with the first data from RHIC p|, 
thus justifying the use of a simple geometric saturation criterion in central collisions. 
To what extent the saturation approach is applicable in non-central collisions needs 
careful further analysis and will not be discussed in this paper. 

In 0, predictions for the measurable final state transverse energy were also presented. 
However, only longitudinal expansion in the plasma phase, using boost-invariant scaling 
hydrodynamics , was considered. Based on a comparison with a systematic study of 
the effects of transverse expansion on top of the longitudinal scaling expansion [T^, it 
was concluded in that a dramatic reduction of transverse energy due to the energy 
loss through the pdV work in the fast initial longitudinal expansion is expected: the 
estimated final dEx/dy at ?/ = is only 1/3 (1/6) of the initially released of the 
minijets at RHIC (LHC). It is clear, however, that the pressure drives matter also 
into the transverse direction. Therefore, we expect that the reduction factors of E^ 
mentioned above should be somewhat smaller when a more realistic hydrodynamical 
treatment is applied. 

The main goals of this paper can be stated as follows: First, we want to improve 
the estimates ^ of the final state particle multiplicities by including a realistic list of 
hadrons and resonance decays in the final state. As a result, each particle will carry 
on the average more than four units of entropy, and the conversion factor between the 
total and charged particle multiplicity will be smaller than 2/3. Second, the previous 
estimates of the measurable final state transverse energies are improved as we are 
now using a hydrodynamic description with transverse expansion included. 

As in PI, we use the saturated minijet system as the initial condition for the hy- 
drodynamic expansion, including now also transverse expansion on top of the scaling 



longitudinal expansion. Our study can be applied in the central rapidity region |TT 
but we cannot study the rapidity dependence of the observables. We focus on central 
or nearly central collisions, for which the final state can be taken to be azimuthally 
symmetric. The observables in non-central collisions, such as elliptic flow require 
a proper treatment of an azimuthally asymmetric expansion and will be discussed 



elsewhere |jT3| 



At collider energies the minijet calculation predicts a QCD plasma in the initial state 
and the hydrodynamic expansion brings the system through mixed phase to a phase 
of hadron gas which finally decouples to noninteracting hadrons. Particle spectra are 
computed by folding the flow at the decoupling with the thermal motion of hadrons 



according to the Cooper- Frye decoupling procedure |]T3[. Spectra of all hadrons and 
hadron resonances up to S(1385) are calculated, followed by the resonance decays to 
obtain the full spectra of stable (against strong interactions) particles. The transverse 
energy and multiplicity densities, either in the pseudorapidity rj or the real particle 
rapidity y, are calculated from the transverse momentum spectra of final particles. 
Predictions for the global observables, the number of (charged) particles and the trans- 
verse energy are made for the energy range from the lowest RHIC energy y/s = 56 
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GeV (Au-Au) to the full LHC energy ^/s = 5500 GeV (Pb-Pb) with and without an 
effective centrality cut of 6% applied. Comparison with the first data from RHIC on 
particle multiplicities is shown. Also comparison with our previous results is made 
both for the multiplicities and for the transverse energies. 



2 Initial conditions 
2.1 Initial energy density 

As discussed in 0, the average saturation scale below which the further parton pro- 
duction is inhibited in a central AA collision is determined from a saturation criterion 
which equates the effective total transverse area of minijets to the effective nuclear 
transverse area, 

Naa{po, v^, 0, AF = 1) X ^ = TiRl. (1) 

Po 

The number of minijets produced above a transverse momentum scale po ^ Aqcd into 
a central rapidity unit AY in an AA collision with an impact parameter b = and 
cms-energy y/s is computed as 

Naa{Po, v^, 0, AY) = 2TAA(0)ajet(po, v^, AF, A), (2) 

where Taa{0) ~ A'^/{7[R\) is the standard nuclear overlap function 0] and a^^t is the 
perturbatively computable minijet cross section with a rapidity acceptance AY. In 
lowest order (LO) 

c^ict{po,Vs,AY,A) = K-Y^ dpl,dyidy2Xifi/A{xi,Q'^)x2fj/A{x2,Q'^)—^ , (3) 

9. 9. 9 

where the minijet rapidities are yi^2- The fractional momenta of the colliding partons 
are xi,2, and the scale is chosen as Q = Pt-, the transverse momentum of produced 
minijets. The parton distributions fi/A{x,Q^) = Rf{x,Q^)fi{x,Q'^) contain nuclear 
effects (shadowing) in Q^) as given by the DGLAP analysis EKS98 |T^. More 



details of the treatment of various subprocesses can be found e.g. in [0. The next- 
to-leading order (NLO) contributions to the number of minijets are simulated by a 
i^'-factor taken from an exact NLO calculation of the minijet transverse energy in 
Ref. [0] at the full RHIC and LHC energies. As the parton distributions in the free 
proton, we use the GRV94 parton distributions [|T7|. It should be noticed that the 
magnitude of the K-iactor depends on y/s and on the parton distributions used [0. 

The saturation scale is obtained as the solution ps^t = Po(\/s, ^) of Eq. (|^). After 
this we compute the transverse energy carried by the produced quanta, 

^T^feat, V^, 0, AY) = TAAiO)a^et{ET)iPs.U v^, ^Y, A), (4) 
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where a^ctiEr) is the first Et moment of the minijet Et distribution. For massless 
partons Et = E sin 6 = pt and the rapidity acceptance AY is implemented by defining 
(here in LO) Et = Pr[0(l/i ^ AF) + 0(?/2 ^ AF)]. The details of the computation of 
(Tjet(-Z?T) in LO can be found in and the extension to the NLO in 0]. 

The determination of Psat also offers a possibility to determine the average initial 
formation time of the QGP as Xj = 1/psat- After this, we can form the average initial 
energy density at z = by applying Bjorken's estimate 0, 

where the rapidity y of the minijet is assumed to be equal to the space-time rapidity 
?/s = I ln[(t + z)/{t — z)] . This implies that the spatial width of a slice corresponding to 
a rapidity interval AY around 2; = is given as Az = TiAY. For hydrodynamics, how- 
ever, we need the energy density with a transverse profile, e = e(s). The simplest way 
of obtaining this is through decomposing T^^(O) = / cP sTa{s)Ta{s) , where the nuclear 
thickness function is obtained from the nuclear density as Ta{s) = J dznAiVs"^ + z^). 
We then arrive at the local energy density profile in the transverse plane 

e(s) = e(s, z = 0) = [T^is)]' ■ ajet(i5T)(Psat, v^, AY, A) ■ (6) 

We use Woods-Saxon nuclear density profile n^(r) with central density no = 0.17 fm~^, 
surface diffuseness d = 0.54 fm and nuclear radius Ra = 1.12A^^^ — 0.86A^^^^. At the 
center, Ta{0) ~ 2i?^no. Notice the azimuthal symmetry as only central collisions are 
considered. 



2.2 Centrality selection 

To have a realistic comparison with the experimental results, we need to simulate the 
centrality selection of the experiments. We do this by considering central collisions of 
an effective nucleus, A^s < A, determined by the centrality selection in the following 
way: We write the total inelastic cross section of the A A collision as cr^'^iy/s) = 
J d'^b[l — exp{—af^ {^/s)TAAi^))]■ Let cr^'^iK) be the inelastic cross section for collisions 
with b < br corresponding to a fraction r (lOOr = n%) of the total inelastic cross section. 
We then have 

r^f/ = ^f/ibr) = d%[l - eM-<TAA{h))] ^ vrfo^ (7) 

The latter approximation is very good for r < 0.1. The average number of participants 
corresponding to the centrality selection of a fraction r can then be estimated as 

{N^.n)r-^ t d%N^^,,{b), (8) 

TTO: JO 
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[GeV] 


^Mcff 


Psat [GeV] 


Ti [fm] 


a{ET) [mbGeV] 


K 


af^ [mb] 


56 


197/177 


0.93/0.92 


0.21/0.22 


40.21/41.67 


2.3 


35 


130 


197/178 


1.08/1.06 


0.18/0.19 


65.17/67.38 


2.3 


39 


200 


197/178 


1.16/1.15 


0.17/0.17 


83.54/86.33 


2.3 


42 


500 


197/178 


1.34/1.32 


0.15/0.15 


131.7/136.0 


2.0 


48 


500 


208/188 


1.35/1.33 


0.15/0.15 


129.5/133.7 


2.0 


48 


1500 


208/189 


1.58/1.56 


0.12/0.13 


217.0/223.5 


1.6 


56 


5500 


208/190 


2.03/2.01 


0.10/0.10 


468.4/481.7 


1.6 


67 



Table 1: The initial conditions for the cms-energies and nuclei considered. The effective nuclei 
Agff correspond to the 6% centrality selection computed with the inelastic cross sections af^ 
as explained in the text. Within a column, the numbers on left correspond to central AA 
collisions and those on right are for central ^efr^cfr collisions. For the highest and lowest 
energies, the i^-factors are taken to be the same as computed in for ^/s = 200 and 5500 
GeV, and for ^/s = 500 GeV in between them. These factors are already included in the 
numbers for {Et)- 

where the number of participants in a collision at an impact parameter b is computed 
from A^part(b) = 2/ (PsTA{\h — s|) [1 — exp{—af^ {^/s)Ta{s))] for symmetric collisions. 

We then find A^q by requiring the number of participants in a central collision of two 
nuclei with mass number Aes to be equal to (A^part)r- Motivated by the first PHOBOS 
results H] we apply here a centrality selection of 6%, which leads to A^s = 177. . . 178 
for Au-Au collisions at RHIC, and to A^s = 188 . . . 190 for Pb-Pb collision at the LHC. 
The information for computing the initial transverse energy profiles is given in Table 1. 

3 Expansion and final particle spectra 

We treat the expansion hydrodynamically. The high initial density with high collision 
rate is one factor supporting the use of hydrodynamics. However, for hydrodynamics 
to apply, the matter should be close to thermal. This is the case with the initial 
conditions which we obtain from the minijet calculation. The system is close to thermal 
equilibrium in the following sense: Using = 1/psat as the time when all the final 
partons have been produced, we can calculate both the initial parton density and 
energy density from the calculated number of partons and transverse energy in 
unit rapidity. For a thermal system these should yield the same initial temperature, 
Tj. Assuming an ideal equation of state of massless partons, we test this by solving 
first Ti from and then calculate the number density from Tj. The result deviates 
from the initial minijet density negligibly at RHIC energies and by few percents at 
the LHC energies, see the curves in Fig. |l[ For the hydrodynamic expansion the 
key quantity is the ratio of pressure to energy density. For massless particles with 
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spherically symmetric momentum distribution this is 1/3 independent of the radial 
shape, the dependence on |p|, of the distribution. Also the main part of the contribution 
to the pressure comes from momenta around the mean value. The tail, where the initial 
deviation from the thermal distribution is probably the largest, contributes less. 

We should also mention that a reasonably good description of the experimental 
data on elliptic flow in terms of hydrodynamics \T2\ gives further confidence in the 
hydrodynamic approach. 

We are here concerned with the expansion of central rapidity region at collider en- 
ergies. Independent of the exact shape of the rapidity distributions they are expected 
to be quite fiat for |?/|<1. The longitudinal expansion should then be well described as 
scaling expansion, Vz = z/t, which means that the flow rapidity equals the space-time 
rapidity Vs = \ ln[(t + z)/{t — z)]. Assuming azimuthal symmetry, as is the case in zero 
impact parameter collisions, reduces the hydrodynamics to 1+1 dimensional problem 
leaving only the radial expansion to be solved numerically [|18] . 



The extra ingredient which is needed for the hydrodynamic calculation, is the equa- 
tion of state (EOS). We use here an EOS A from Ref. |T^, which was employed in 
a hydrodynamic study of hadron spectra at SPS energy. The quark-gluon-plasma 
phase is treated as an ideal gas of massless gluons and two flavours of quarks and 
antiquarks. In hadronic phase hadrons and hadron resonances up to S(1385) are in- 
cluded and the repulsion among hadrons is described through a mean field parameter 
K = 450 fm^MeV. At the phase transition temperature, Tc = 167 MeV for zero 
net baryon number, matter with densities between those of plasma and hadron gas is 
assumed to be in an equilibrium mixed phase. 

The spectra of all hadrons and hadron resonances up to S(1385) are calculated when 
the density drops to values where the mean free paths become similar to the dimension 
of the system. The decoupling condition is defined in terms of the energy density which 
essentially fixes the temperature since the values of baryon chemical potential are small 
at collider energies in the central rapidity region. At the decoupling we have T^ec — 120 
MeV. It can be expected that the hadrons with smaller cross sections decouple at 
earlier time. Similarly, the flavour changing cross sections are small and are known to 
decouple earlier. However, for the transversally integrated quantities considered here, 
these effects are small and we have taken the same decoupling condition for all particles 
and assumed both kinetic and chemical equilibrium at the decoupling. For the folding 
of flow and thermal motion we have employed the Cooper and Fry |jl^ prescription to 
calculate the spectra. 

To obtain the final stable (against strong interactions) hadron spectra, two and three 
body decays of all hadron resonances pO|] have been included. The details of the particle 
spectra, like their dependence on the mass number and energy of colliding nuclei, will 
be considered in a separate publication . 
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4 Multiplicities and transverse energy 



From the treatment of the decouphng and resonance decays, we obtain particle spec- 
tra dN/ d'^pTdy which are boost invariant and azimuthally symmetric by construction. 
These symmetries are utihzed in the definitions of multiphcities and transverse energies 
below. 

1. The rapidity distribution of the total number of particles, the total multiplicity. 



dy 



dN, 



y=0 



(9) 



dpTdy 

where the index i runs through all particle species. Due to the boost invariance dN/ dy 
is independent of y and equals the rapidity distributions averaged over a rapidity bin 
Ay. 

2. The pseudorapidity distribution of the total multiplicity is defined as 



dN 



where y 



dprdr] 
arsinh(^^ sinh?7), and 



dNi{pT,y) 
dprdy 



J{v,pt) 



dy_ 



p_ 



Pt 



(10) 



:iii 



Ti 



with the transverse mass m^- 



p"^ + mf. The averaged total multiplicity of particles 
in a pseudorapidity bin At] symmetric around t] = can then be defined as 



dN 
drj 



1 



dN 



A»? A?7 jaj7 d7] Arj 



dN 



dpxdy 



prp Ari 

arsinh( sinh ), (12) 

y=Q rriTi 2 



where in the last step the boost invariance was used. To obtain the charged particle 
multiplicity, one simply excludes all neutral particles from the sum. When applying 
these results here, we will choose Ar] as |?7| < 1 in accordance with the PHOBOS 
experiment 

3. We will be interested in two transverse energy quantities. First, in direct corre- 
spondence with the experimental definition of the Et as the energy in a calorimeter cell 
at certain rj or scattering angle 6, we define Et = E sin 6, and get the Et distribution 
from the calculated particle spectra as 



dET 
drj 



= dpTY. 



dpTdrj 



-E' 



Ti 



dprY. 



dN 



dpTdy 



y=0 



Pt- 



(13) 



Note that again due to the boost invariance the result is independent of rj, so averaging 
over a pseudorapidity bin Arj gives conveniently 



dEj 



drj 



dET dET 
Art A?7 Jat] dr] dr] 



(14) 
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The other transverse energy distribution we wish to study is defined as 



dET 



dNi 



(15) 



dy y=o 



Y dprdy y=o 



rriTu 



where we have used Et = rriT as is the case at ?/ = 0. 

5 Results 

The total particle multiplicity dN/dy at ?/ = 0, computed from Eq. (j^) by using the 
particle spectra, is shown in Fig. |l]. We consider Au-Au collisions (triangles) at RHIC 
energies ^/s = 56, 130 and 200 AGeV and, as an interpolation between RHIC and 
LHC, also 500 AGeV. Pb-Pb collisions (squares) are considered at the LHC design 
energy ^/s = 5500 (LHC) and also at 1500 and 500 AGeV. The open symbols show 
the results for central collisions at b = 0. The filled symbols stand for the results with 
a 6% centrahty selection, described as central collisions of effective nuclei A^s < A, 
see Table 1. As seen in the figure, the multiplicities are reduced by about 10% by this 
centrality cut. 

The multiplicities for central collisions from our previous work are shown as 
lines and they should be compared with the open symbols. In [^, the final state 
multiplicity was obtained either by converting the initial state parton multiplicity into 
entropy, leading to Nf{Ni) = Si{Ni)/A = (3.6/4) * 1.383^0-922(^^)0.383 (g^jj^ jj^gg)^ 
or by converting the initial state energy density into entropy, leading to Nf{ETi) = 
Si{ETi)/4: = 1.16y4°-92(A/i)°-'^° (dotted lines). As observed in @, these two results are 
close as the average initial energy per particle is very near thermal. 

The multiplicities obtained for central collisions in the present study tend to be 
above EKRT at RHIC energies but below at the LHC energies. As mentioned in the 
introduction, this is due to different but partly compensating effects that appear first 
in the conversion of the initial energy density into the initial entropy and then the 
entropy into the final state multiplicity: 

• In Pi , only the gluonic degrees of freedom were used to obtain temperature 
but here we choose to include also two fiavours of quarks and antiquarks. The 
resulting increase in the initial entropy density (computed from e) enters through 
s ~ g^^'^e^^'^, with g as the effective number of the degrees of freedom, causing a 
relative increase of about 23%. 

• In comparison with 0, where no transverse profiles were included, the total 
rapidity density of initial entropy, dSi/dyg, is reduced by about 5% due to the 
local conversion e —>■ s. 

• During the hydrodynamic evolution stage, some entropy is generated at the phase 
transition and by the numerics. These effects are quite small, e.g. at y/s = 500 
AGeV the increase in entropy is about 3%. 
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Figure 1: Total particle multiplicity dN/dy for ^ = 56, 130, 500, 1500, 5500 ^GeV. The open 
circles are for central Au-Au collisions, the open squares for central Pb-Pb collisions. The 
filled symbols are for central collisions of jdefjAgfj collisions, corresponding to a 6% centrality 
cut (see Table 1). The results for central Au-Au and Pb-Pb collisions from Q], computed 
either from the initial multiplicity (solid), or from the initial energy density (dotted) are 
shown respectively. 

• In a realistic hadron resonance gas of massive particles at freeze-out (now T^cc = 
120 MeV), the average entropy per particle is larger than four, Sf/Nf ^ 4.68 
with the list of hadrons used. This decreases the multiplicity from by about 
15 %. 

Combining the effects above we get an estimate of dN/dy in terms of Nf{Ej'i) of 
EKRT as dN/dy = 0.95 * 1.23 * 1.03 * (4/4.68) * NfiEn) = 1.03 * NfiEn). This 
explains the small difference between the dotted lines and open symbols in Fig. Q at 
^/s = 500 y4GeV, where the i^'-factor applied coincides with the constant i^-factor of 
EKRT. Larger differences at other ^/s result mainly from the different i^'-factors, listed 
in Table 1. 

The corresponding charged particle multiplicities, averaged over two central units of 
pseudorapidity, dN/dr]\\r,\<i, computed from Eq. ([T^), are shown in Fig. |^. Again, the 
open symbols denote the central collisions and the closed ones have the 6% centrality 
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cut included. Otherwise the notation is identical to that in the previous figure. The 
PHOBOS data at y/s = 56 and 130 AGeV |Q] with the 6% centrality selection, is shown 
by open circles. The results (filled triangles) are seen to agree remarkably well with the 
data, given that no fitting was done. We should mention, though, that the PHENIX 
data at y/s = 130 AGeV (not shown in the figure) lies on the upper edge of the 
corresponding PHOBOS point. As indicated by the discussion above, the uncertainties 
in the theoretical computation are larger than the error bars of the data. 
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Figure 2: Charged particle multiplicity (iA^/(i?7|[^|<i averaged over a pseudorapidity bin —1 < 
r/ < 1, computed from Eq. (^). The cms-energies, the nuclei, the centrahty cuts and the 
symbols are as in Fig. |l[ The PHOBOS data (open circles) are shown with systematic error 
bars, the small statistical error bars are not visible. The solid and dotted lines are again the 
predictions presented in for the central collisions with the same notation as in Fig. |l[ 

Again, the EKRT results for dNch/dr] from 0] are shown for the central collisions 
(b = 0) by the solid and dotted lines, scaled down from those in Fig. ^ by a factor 
2/3 to account for the conversion to charged particles and a factor 0.9 to account 
for the conversion from particle rapidity to pseudorapidity and averaging over the 
pseudorapidity bin Ir^l < 1. The slight decrease in the ratio of the present results (open 
symbols) to those of EKRT as compared with the same ratio for total multiplicities, see 
Fig. m, is mainly due to the inclusion of resonance decays, which reduces the effective 



10 



ratio A^chMot from 2/3 to 0.6. 

Note that if the EKRT results are multiphed by a further factor 0.9 to account for the 
6% centrahty selection (or if A^s is used [^), they agree very well with the first RHIC 
data from PHOBOS. The effects discussed above, related to the computation of the 
initial entropy and conversion of the final state entropy to the final state multiplicity, are 
to a good approximation constant multiplicative factors at all energies. These factors 
and the varying iT-factors from the pQCD computation can be easily included in the 
EKRT saturation model which thus provides an effective approach for the computation 
of the final state multiplicities in central collisions. If, however, particle spectra need 
to be computed, a more detailed approach, such as the present one, is necessary. 
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Figure 3: The total transverse energy distribution dEx/dy at y = computed from Eq. ( p!5| ) 
for Au-Au and Pb-Pb central collisions at various energies with and without a 6% centrality 
cut. The crosses show the initial E^"^ in Ay from minijets at the saturation scale. The solid 
lines are the prediction of E^^ from and the dotted lines are the initial E^^ of § for 
central collisions. Otherwize the notation is the same as in Figs. 1 and 2. 



In the computation of the final state transverse energies the transverse expansion 
effects play an important role. Next, we plot the final state dErp/dy at ?/ = from Eq. 
(p!5|) in Fig. |^ for the same energies, nuclei and centralities as in the previous figures. 
To get an idea of how much of the initial transverse motion of minijets is converted 
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into longitudinal motion in thermalization and, especially, by the pdV work in the 
expansion, we show E^^ of minijets at saturation by the dashed line, computed from 
Eq. (H). The EKRT results for the final state transverse energy computed in as 
Erpf = 7TR\e{Tf)Tf are shown by the solid lines for central A A collisions. In |^ the time 
Tf, at which the evolution was terminated, is the time of reaching the energy density 
ec = e(rj) estimated from the results of The dashed lines show the initial of 
1^ in central collisions with a constant factor K = 2 included at all energies. Notice how 
the reduction of transverse energy between the initial and final state becomes smaller 
with the inclusion of the transverse expansion in the hydrodynamical description. As 
seen in the figure, the initially released Et is then reduced by a factor 2.7 at the full 
RHIC energy, and 3.6 at the LHC 

As the last item here, we study the measurable total transverse energy pseudorapidity 
distributions of the final state hadrons. Fig. ^ shows our predictions for dEx / drj\\j^\<i 
averaged over the pseudorapidity bin \r]\ < 1, computed from Eq. (|1^) for the same 
energies, nuclei and centralities as before. 

We have also checked the effect of the decoupling temperature T^ec on the multiplic- 
ities and transverse energies. If Tdec is changed by ±20 MeV, the transverse energy 
changes respectively by ±2 . . . 3% both at RHIC and at LHC energies. The correspond- 
ing changes in the total multiplicities are ±1% both at RHIC and at LHC. 



6 Discussion and conclusions 

We have made predictions for the measurable multiplicities and transverse energies in 
central and nearly central Au-Au and Pb-Pb collisions at collider energies. Hydro- 
dynamics with transverse expansion and a realistic equation of state, supplemented 
by resonance decays in the final state has been applied. Contrary to the situation 
at the SPS, where the measured hadron spectra were needed to constrain the initial 



conditions [^], we now compute the initial conditions from perturbative QCD, supple- 
mented with the saturation of produced partons 0. This is an essential improvement. 
For the final state multiplicities the results agree very well with the first PHOBOS 
results. It should be emphasized that no fitting or fine tuning was done in order to 
arrive at this result. 

The multiplicities obtained in this study are quite close to those of but the final 
state transverse energies are now about 20. . . 40 % larger than those predicted in In 
comparison with EKRT, there are different compensating effects but the difference is 
mainly due to the inclusion of the transverse flow which reduces the transverse energy 
loss but does not change the entropy. 

It is difficult to estimate the uncertainty in the saturation approach since the satu- 
ration criterion itself contains an undetermined constant of order 1, containing group 
theory factors and powers of and taken equal to unity in H as well as here. This 
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Figure 4: Prediction for the total transverse energy distribution dEx/dr] at = as given 
by Eq. (|T3| ) for Au-Au and Pb-Pb central collisions at various energies with and without 
a 6% centrality cut. The notation is the same as in the previous figures. The curves are 
parametrizations to guide the eye. 



uncertainty affects the absolute normalization of the results, but we note that data on 
multiplicity at one value of ^/s and A is sufficient to fix the constant, therefore allowing 
predictions for absolute values of the global quantities at different ^/s or A. The real 
predictions of the saturation model are rather the scaling laws as a function of a/s and 
A for the global quantities in central or nearly central collisions. The predicted ^/s 
scaling seems indeed to be compatible with the first PHOBOS data 0. 

To improve these calculations further one should use the local saturation criterion ||2^ 
to obtain the initial profiles. The approach of as such, however, is not directly 
applicable in this context since there one considered only the saturated part of the 
initial conditions. For the hydrodynamics one needs also to take into account the 
region beyond this, i.e. to consider the tails of the matter distributions. The transition 
from saturated to non-saturated initial conditions is a delicate issue concerning the 
interplay between geometry and dynamics of the collision and to be studied in the 
future. 

As noted above, the y/s scaling of multiplicity in the saturation model agrees well 
with the results of PHOBOS. This gives more confidence in extrapolating the calcu- 
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lations towards LHC energies. Our prediction for the average charged particle mul- 
tiphcity in central Pb-Pb collisions at y/s = 5.5 TeV is about 2560. This is slightly 
smaller than the prediction of 0, as the energy dependence of the ii'-factor is taken 
into account. The charged particle multiplicity is thus clearly less than the LHC design 
value 8000. 

While our results for the multiplicity agree with the PHOBOS data, they are a 



little below the corresponding PHENIX data [|^. For the Et our results are a little 
above the preliminary PHENIX data In this study we have not tried to tune 

the calculations in order to arrive at some desired, measured, numbers on either N or 
Et- Rather, we have shown that the initial conditions from pQCD minijet calculation 
featuring a dynamical saturation scale and followed by hydrodynamic evolution with 
transverse expansion, lead to results which are very similar to the measured values. 
More precise tests, e.g. in terms of the transverse momentum spectra of different 
particles, are needed to learn more details of the transverse dynamics and to be able 
to pin down the properties of the initially produced matter. 
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